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Abstract — There  is  a  consensus  that  cooperative 
communication  is  a  key  technology  for  the  upcoming 
generations  of  mobile  communications.  Whether  under 
the  Coordinated  Multipoint  (CoMP),  Device -to-Device 
(D2D)  communication  or  Internet  of  Things  (loT)  scopes, 
there  is  little  doubt  that  a  more  decentralized  network, 
where  several  communicating  agents  smartly  cooperate 
to  improve  data  rate  and  reliability,  is  fundamental  to 
cope  with  the  always  scarce  resources  of  energy  and 
bandwidth.  When  relaying  stations  are  deployed  to 
increase  signal  coverage  and  to  provide  cooperative 
diversity,  the  concatenation  of  two  or  more  point-to-point 
MIMO  channels  impose  a  restriction  on  conventional 
algorithms  for  channel  estimation.  Although  single-relay, 
two-hop  Amplify-and-Forward  (AF)  systems  have  been 
reasonably  investigated,  there  are  enough  evidences  that 
more  transmission  hops  will  be  a  common  scenario  in  the 
next  years.  This  work  proposes  a  semi-blind  receiver  for 
the  task  of  joint  symbol  and  channel  estimation  in  a  three - 
hop  AF  MIMO  system.  Resorting  to  a  Nested  PARAFAC 
tensor  model,  new  estimating  equations  were  derived, 
along  studies  on  the  uniqueness  and  i dent ifi ability 
conditions  concerning  the  aforementioned  receiver. 
Simulations  corroborate  the  validity  of  the  receiver  as  an 
effective  option. 

Keywords — Amplify-and-Forward  ,  Cooperative 
communications,  Relaying,  Tensor  decomposition. 

I.  INTRODUCTION 

hi  wireless  communications,  relay  stations 
havebeen  extensively  deployed  to  increase  the  signal 
coverage  in  a  broadcast  transmission  [15].  hi  recent  years, 
relaying  also  has  been  used  to  provide  digital  receivers 
with  what  is  called  cooperative  diversity,  where  a  relay 
might  work  with  single -antenna  transmitters  to  emulate  a 
virtual  array  of  antennas.  In  this  sense,  the  benefits  of 
spatial  transmit  diversity  could  be  obtained  without  a 
increase  of  the  number  of  transmit  antennas. 

The  simplest  of  the  relaying  protocols  is  called 
Amplify-and-Forward  (AF)  As  its  name  suggests,  the 
incoming  signals  are  amplified,  and  then  forwarded  to  the 
intended  destination.  The  main  advantage  of  this  non- 


regenerative  process  is  that  no  complex  decodingis 
performed  at  the  relay,  so  its  hardware  and  software 
complexities  arereduced  with  comparison  to  the  so-called 
regenerative  protocols,  such  as  the  Decode-and-Forward 
(DF),  favoring  its  mass  implementation  [15-16].  As  a 
drawback,  due  to  its  simple  operation,  noise  and 
interference  are  also  amplified,  limiting  the  benefits  of 
relaying. 

Another  disadvantage  of  the  AF  protocol  is 
related  to  the  task  of  channel  estimation.  Besides  the 
usual  need  of  the  Channel  State  Information  (CSI)  for 
symbol  detection,  many  system  optimization  techniques 
demand  the  knowledge  of  the  various  channels  that 
compose  a  relaying  network  [1,  8-10].  The  absence  of  a 
powerful  processing  unit  at  the  relay  leaves  all  this  work 
to  the  destination  node,  which  is  usually  inapt  to 
dissociate  the  cascaded  source-relay  and  relay -destination 
channels  under  conventional  transmission  protocols. 

The  issue  of  channel  estimation  in  two -hop  AF 
MIMO  systems  have  been  addressed  in  few  state-of-the- 
art  works  [5-7,14]  by  using  training  sequences  and  by 
[11,  12]  in  a  semi-blind  fashion.  However,  when  it  comes 
to  the  multi-hop  systems,  with  more  than  one  serial  relay 
in  the  communication  link,  there  is  usually  a  lack  of 
related  works.  This  is  particularly  concerning  since  the 
future  of  wireless  communications  points  towards  ever 
more  decentralized  networks,  such  as  Device-to-Device 
(D2D)  communications  technologies,  where  longer  chains 
of  connected  relaying  devices  can  be  expected.  For  a 
three-hop  (two  relays)  link,  [2]  have  resorted  to  the 
combination  of  PARAFAC  and  Tucker  tensor  models  to 
estimate  the  channel.  However,  [2]  resorts  once  again  to 
the  use  of  pilot  symbols,  disregarding  the  benefits  of  the 
joint  estimation  of  symbols  and  channels. In  [18]  the 
authors  generalized  the  works  of  [11,12]  for  any  number 
of  hops. 

This  work  proposes  a  three -hop  AF  MEMO 
system,  where  both  relays  and  the  source  node  applies  a 
Khatri-Rao  Space-Time  (KRST)  coding  over  the  signals 
prior  to  their  transmission.  In  this  way,  the  receiver  can 
arrange  the  signals  into  a  5th -order  tensor  following  a 
Nested  PARAFAC  model,  from  which  the  three-hop 
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channels  and  symbols  can  be  estimated  by  a  closed  -form 
algorithm.  Notably,  this  proposed  tensor  model  is  an 
expansion  of  the  4th-order  model  ,holding  similar 
properties/advantages.  Studies  on  the  uniqueness  of  the 
proposed  model  and  and  on  identifiability  conditions  for 
the  proposed  receiver  are  also  carried  out,  along  the 
presentation  of  its  computational  complexity  and 
performance  in  terms  of  Bit  Error  Rate  (BER)  and 
channel  Normalized  Mean  Square  Error  (NMSE). 

The  main  contributions  are: 

1.  Derivation  of  new  unfolded  equations  for 
the  three-hop  system,  for  both  signals  and  noise  terms; 

2.  Proposal  of  a  theorem  on  the  essential 
uniqueness  of  the  systemmodel; 

3.  Discussion  on  the  computational  complexity 
and  on  the  achievable  spectral  efficiency  of  the  proposed 
method. 

1.1  Notations  and  fundamentals  of  tensors 

Scalars,  column  vectors,  matrices,  and  tensors 
are  denoted  by  lower-case  (x),  boldface  lower-case  (x), 
boldface  capital  (X),  and  calligraphic  (X)  letters, 
respectively.  XT,  X*,  X*,  X;.  and  X.m  are  the  transpose, 
the  conjugate,  the  pseudoinverse, the  Ith  row,  and  the  mth 
column  of  X  £  CLxM,  respectively.  Dn(X)  stands  for  the 
diagonal  matrix  formed  from  the  elements  of  Xn..  Given  a 
third-order  tensor  X  £  C,XJxK ,  with  entry  xijk,  the 
matrices  XJKxI,  XKIxJ  and  XlIxK  denote  tall  mode-1, 
mode-2  and  mode-3  unfoldings,  with  xijk  = 

^JKx,\k-l)J+j,i  KIXj\i-i)K+k,j  )/+ik' 

The  vec(.)  and  unvec(.)  operators  are  defined  by 
x]KI  =  vec(xJKxI)  £  CJK,X1 
~  XJKxl  =  unvec(x;ffl). 

1.2  PARAFAC  and  Nested  PARAFAC 

decompositions 

A  PARAFAC  decomposition  [4,13]  of  a  third- 
order  tensor  X  £  C,x]xK,  with  rank-/?  and  matrix  factors 
(A,  B,  C),  will  be  noted  lA,  B,  C;  /?]  .  Tall  and  flat  mode-1 
matrix  unfoldings  of  X  are  respectively  given  by 

X]KX,  =(C«B)Ar  =(XIXJK)T, 
where  o  denotes  the  Khatri-Rao  product.  Similar  mode-2 
and  mode-3  unfoldings  can  be  obtained  by  permuting  the 
factor  matrices. 

A  4th  Nested  PARAFAC  model  [17,18]  can  be 
seen  as  a  generalized  model  w.r.t.the  PARAFAC.  hi  this 
model  a  4th-ordei  tensor  can  be  reduced  to  two3rd-order 
tensors  that  follows  a  PARAFAC  model,  by  associating 
two  of  their  original  modes  into  a  single  one.  The  matrix 
factors  related  to  the  associated  modes  are 
X  =  lA,  B,  C;  Rj, 


C  =  lD,E,F;/?2l. 

where  C  can  be  any  matrix  unfolding  of  C.Note  that  C  in 
(2)  is  any  unfolded  form  of  C  represented  in  (3). 


II.  THREE-HOP  SYSTEM 

A  one-way  three -hop  relay  system  is  depicted  in 
Fig.  1.  The  source  (5)  node  wants  to  transnit  its  signal  to 
the  destination  (D)  node  with  the  aid  of  unidirectional 
relays  and  /?2.  The  communication  channels  H^1-1  £ 

£M1XMS  ^  £  £M2XMt  an(J  pj(fl2D)  £  £MDXM2  are 

considered  to  flat-fading  and  time -invariant  during  a 
transmission  block 


Ms 


M, 


M2 


Md 


Fig.l:  Three-hop  one-way  communication 


Prior  to  the  transmisstion,  the  source  node 
applies  a  Khatri-Rao  Space-Time  (KRST)  coding  on  its 
symbol  matrix  S  £  £NxMs  using  the  code  matrix  cPxMs . 
Just  like  in  ,  N  is  the  data-stream  (number  of  symbol 
vectors)  and  P  is  the  spreading  source  code  length.  Thus, 
in  the  first  hop  of  Fig.l,  the  signals  transmitted  from  S'  to 
/?!  become 

x'im1xpn  n  w  ■ 

In  the  second  hop,  R1  applies  the  Khatri-Rao  ST  coding 
over  /,  and  R±  transmits  its  signal  to  R2.  The  signals 
arriving  at  this  second  relay  are 

W(/?z)  =  H|Ti,,2)('w(Tri  tG(#i)f 
vvm2xjpn  n  VvvPJVxM1  u  J  ■ 

Equation  (5)  is  intentionally  written  in  the  form 

of  (4),  such  that  there  is  a  correspondence  between  the 

matrices  that  compose  both  signal  models.  For  one,  G^1^1 

has  equivalent  role  of  C  on  the  R/s  side,  while  WpJy^Ml 

and  S  represent  the  signals  prior  to  the  ST  coding  at  the 

first  relay  and  at  source,  respectively.  The  entris  of 

^m^xjpn  could  be  stored  in  the  4th  -order  tens  or  £ 

£M2xjxpxn  such  as 

u,(R2)  =  VMl  h(RiR2)  n(Ri)yMs  h(sRi)  (1) 

wm2  ,j,p,n  ^-‘m1  =  1  Lm2  ,m1  &j  ,m1  =  l  ILm1 ,171$  '-'p.ms^n.ms 

that  follows  a  Nested  PARAFAC  decomposition 
However,  without  needing  to  store  the  data,  R2  forwards 
the  signals  to  the  destination  node  (last  hop  of  Fig.  1). 

Before  the  forwarding  process  by  R2,  this  relay 
applies  a  new  orthogonal  KRST  coding  (similarly  to  (4) 
and  (5))  with  a  code  length  K,  such  that  after  the 
transmission  through  H^2®^  we  have 

X  =  H  (R2D)  fw  <’7?2-1  oG^2^7 

AMDXKJPN  n  VVvypjVxM2  u  >  ■ 

(2) 
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Due  to  the  third  hop  (third  KRST  coding),  with 
the  addition  of  the  dimension  of  length  K,  the  signals 
arriving  at  the  destination  node  can  be  arranged  in  a  5th- 
order  tensor  X  G  £MDxkxjxpxn  the  scalar  form 

of  this  tensor  can  be  given  by 
m2 

r  =  V  h(R2D)  n(R2)u,(R2) 

mD,k,j,p,n  nmD,m29k,m2  m2,j,p,n 

m2  =  1 

=  yM2  (R2) 

^~>m2  =  1  ,m2y k ,m2 

yMi  h  (RtR2)  (i?i) 

=  1  'Lm2,m1  &j  ,m1 

x  y Ms  h(SR i) 

2-‘ms  =  1  m-L.ms  V.™-s  n,ms- 

Contrasting  (8)  with  (6),  the  nesting  nature  of  the 
Nested  PARAFAC  is  kept  at  a  higher-order  due  to  the 
sequential  KRST  coding  at  each  relay  station.  For  the  4th 
order  model,  the  data  can  be  processed  directly  from  the 
4th -way  tensor  or  through  the  two  nested  3th-order 
tensors,  being  the  latter  the  best  in  terms  of  computational 
burden.  Among  different  forms  of  processing  the 
(received)  5th-order  tensor  described  (8),  the  data  can  be 
arranged  to  conform  to  three  PARAFAC  models, 
X(z) £mdrJ  xpxn ,  21  ( Rz )  e  £MDkx]xms  an(j  2,(r l)  £ 

£MDXKXM1 

^(Z)  =  IzS/xms-C,S;Ms], 

^(*2)  =  [zmSxMi.g(*i).(h(m0)i’;m11. 

Z(Rl)  =  |H(i?2D),G(,?2),(H(/?1,?2)):r;M2|, 

where  X  ^  is  obtained  through  the  order  reduction  of  X 
by  associating  its  three  first  indices  .  The  predilection  for 
working  with  3  these  3rd -order  tensors  rather  than 
directly  with5f,i -order  tensorcomes  from  the  fact  that 
established  tools  and  theorems  dedicated  to  third-order 
PARAFAC  models  are  well-established  and  could  be 
more  easily  applied. 


2.1  Noise  terms 

The  additive  noise  terms  on  relays  R1  and  R2  can 
be  respectively  given  by 


\/(R  D  —  of.th))1, 

vm2xjpn  11  \vm2xpn  u  J  ’ 

y(RD)  =  u(R2D) o  G^2^1T 

These  forms  follow  the  same  unfolding  order  of 
(5)  and  (7).  The  matrices  V^p„,  yjg/pw  and  VmdDxkjpn 
are  unfolded  forms  of  the  noise  tensors  G 

£M1XPXN  y(R2)  e  £M2X]XPXN y(D)  £  £MDxKxJxPxN 

that  affect  the  inputs  of  R1?  R2  and  of  the  node  D, 
respectively. 

Note  that  V^2-1  and  ’V^>^KjPN  also  follow  Nested 
PARAFAC  decompositions,  and  thus,  its  matrix  factors, 
such  as  the  channel  and  code  matrices,  night  be  uniquely 


identified.  That  is  a  interesting  approach,  as  it  would  not 
require  symbol  transnission  for  the  task  of  CSI  retrieval. 


HI.  TRIPLE  KRF  (TKRF)  SEMI-BLIND  RECEIVER 

The  objective  of  the  semi-blind  receiver  installed 
at  the  destination  is  not  only  to  decode  the  symbols 
transmitted  by  the  source,  but  also  to  jointly  estimate  all 
channels  of  the  3-hop  link.  The  use  of  the  4th-order 
Nested  PARAFAC  decomposition  employed  by  [?,  12] 
enabled  the  development  of  3  semi-blind  receivers  in 
total,  although  others  could  have  been  presented  for  the 
same  mdodel.  Perhaps  the  most  interesting  $$  those 
receivers  is  the  one  called  DKRF,  composed  by  two  non¬ 
iterative,  sequential  Khatri-Rao  Factorizations  (KRF), 
which  presented  the  same  performance  as  the  others,  but 
with  less  computational  complexity  .  The  same  idea  can 
be  extended  from  the  two-hop  systemto  the  three-hop  one 
by  adding  up  a  third  KRF  factorization. 

Let  one  write  the  mode-2  unfoldings  of  the 
tensors  X(z),  and  2T(Rl)  respectively  as 


*NMdKJxP 


(z 


(r2) 

Md  KJxM$ 


s)cT, 


<siDKxj  =  (z£U  o  (H<“J)r)(G<*Or.  (13) 

z  %iDXK  =  (h^d>  o  (9) 

The  basic  KRF  factorization  consists  of 
reshaping  each  column  of  a  Khatri-Rao  prod*Uc?  into  a 
matrix,  and  then  approximating  such  rank-one( lrhhtrix  as 
the  outer  product  of  two  column  vectors  .  The  KRF 
algorithm  is  in  Alg.  2 

Here  we  admit  the  principle  that  the  estimation 
of  the  symbols  is  indispensable,  while  the  estimation  of 
the  channel  is  optional.  This  makes  the  TKRF  receiver 
more  flexible,  since  some  of  the  steps  for  the  estimation 
of  individual  channels  can  be  neglected,  resulting  in  the 
reduction  of  the  global  complexity  and  a  relaxation  of  the 
identifiability  conditions.  The  three  (sequential)  KRF  of 
the  TKRF  receiver  are  called  KRF-1,  KRF-2  and  KRF-3, 
and  they  basically  differ  in  the  inputs  and  outputs.  The 
fluxogram  of  the  TKRF  receiver  is  shown  in  the  Fig.  2, 
and  the  description  of  its  inputs  and  outputs  are  given  by: 

•  KRF-1:  Symbol  estimation  (mandatory) 
o  Inputs:  XNMdK]xp  and  C; 

o  Outputs:  ( ZmdIjxms’S) 


•  KRF-2  (optional) 

o  Inputs:  (Z[Rs2lDKXJ,G^y, 
o  Outputs:  (Z^XAfi ,  H (SRl)) 

•  KRF-3  (optional) 

o  Inputs:  (Z^DXjf,G(R2)); 
o  Outputs:  (H^R2D\H^R2^) 
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3.1  Identifiability  and  Uniqueness  Conditions 

In  this  section,  identifiability  and  uniqueness 
conditions  are  derived  to  the  case  where  all  matrices 
(symbols  and  channels)  are  to  be  estimated. 

Theorem  1  (Identifibiality)  Let  the  source  code 
matrix  C  and  the  relay  code  matrices  and  G(■,?2•, 

have  full-rank.  Necessary  and  sufficient  identifiability 
condition  to  jointly  estimate  symbol  (S)  and  channel 
(H(sr'\  H{r^  and  H^r zS))  matrices  are  that 

P  >  Ms,  ]  >M1  and  K  >  M2. 

Proof.  From  the  KRF  algoritm  (Alg.  2(c)@),  the 
necessary  (and  only)  condition  to  apply  the  column 
factorization  is  that  GT(GT)t  =  Ifl.  For  the  existence  of 
the  right  inverse  of  Gr,  if  G  is  full-rank,  then  it  is 
necessary  that  G  be  also  full  column -rank.  Once  the 
TKRF  receiver  calls  the  KRF  algorithm  three  times  (i.e. 
KRF-1,  KRF-2  and  KRF-3),  using  one  of  three  different 
correspondences  at  a  time,  i.e.  G  »  (c,  G^^G  (^?2-1),  then 
the  necessary  and  sufficient  identifiability  condition, 
given  that  all  code  matrices  are  full-rank  by  Theorem  1,  is 
simply  (15). 

Theorem  2  (Uniqueness)  Assume  the  following 
hypotheses:  (a)  the  entries  of  the  channel  matrices  are 
drawn  from  stochastic  processes  with  continuous 
Gaussian  distribution;  (b)  the  identifiability  hypotheses 
and  conditions  from  Theorem  1  are  met;  (c)  Symbol 
matrix  is  composed  of  random  symbols,  and  the  number 
of  data-stream  N  is  much  greater  than  the  number  of 
source  antennas  Ms.  Under  such  hypotheses,  a  sufficient 
but  not  necessary  condition  to  ensure  with  very  high 
probability  the  uniqueness  solution  at  the  output  of  the 
TRKF  receiver  is 

min(MD,Ms,M1,M2 )  >  |M2  -  AfJ  +  2. 

Proof.  Model  equations  (10)  and  (11)  belongs  to 
a  4th-order  Nested  PARAFAC  decomposition  of  the  data 
fully  stored  in  Z^R2\  so  by  using  Theorem  2  from  [12] 
with  due  correspondences,  it  is  possible  to  ensure  the 
essential  uniqueness  of  matrix  factors  in  (10)  and  (11)  if 
the  entries  of  H(Rl,,?2-)  are  drawn  from  a  continous 


Gaussian  distribution.  Once  this  is  the  case,  as  stated  by 
the  hypotheses  of  Theorem  2  in  this  paper,  then  from  [12] 
the  uniqueness  of  the  estimates  of  H^SRl\  H  6*1*2)  and 
H(fl2s)  are  guaranteed  if 

+  fc(H(-s«i))T  -  max(2Mi  -  M2,M2)  +  2, 
feG(R 2)  +  kH(R2D )  >  max(2M2  —  M^Mf)  +  2, 
where  kA  is  the  Kruskal  rank  of  A  . 

Once  G^  and  G <R2>  are  full  column-rank  to 
satisfy  the  identifiability  condition  of  Theorem  1,  then 
feG( Ri)  =  Afy  and  fcG(f?2)  =  M2.  Given  the  hypothesis  on 
the  channels,  we  have  that  k^SRi^r  =  min(Ms,  Mf)  and 

kH;R2D )  =  min (MD,M2).  Eqs.  (17)  and  (18)  respectively 
become 

min (MS,M1)  >  max(M1  —  M2,M2  —  Mf)  +  2, 
min  (Md,  M2)  >  max(M2  —  Mv  M1  —  M2j  +2. 

Finally,  combining  (19)  and  (20)  leads  to  the  condition 
(16). 

Proven  the  uniqueness  of  the  fourth-order  model 
comprising  the  two  three -order  models  (10)  and  (11), 
what  is  left  is  the  need  to  find  the  uniqueness  condition 
regarding  (9).  This  third-order  tensor  model  has  the 
following  Kruskal  condition:  (15) 

k  ;r2)  +  kG  +  k$  >  2MS  T  2. 

zmdkj*ms 

From  hypotheses  (b)  and  (c)  from  Theorem  2,  C 
and  S  (with  very  high  probability)  have  full  column -rank. 
Thus,  kc  =  Ms  and  ks  =  Ms ,  reducing  (21)  to 
k  <r2)  >  2.  This  is  always  tme  if  Z^'L  xM  does  not 

ZM[)KJxMs  D  S 

have  any  zero  column  or  any  pair  of  linearly  dependent 
columns  .  Since  channel  matrices  are  randomly  drawn 
from  Gaussian  distributions,  and  GlRll  and  G <R2>  are  full- 
rank,  then  this  condition  is  satisfied  with  probability  close 
to  one. 

3.2  Computational  complexity 

The  KRF  algorithm  consists  of  multiples  rank- 
one  approximations  using  SVD.  Apart  from  other 
eventual  operations,  such  as  buffering  and  memory 
allocation,  one  may  count  those  SVD’s  as  the  cost- 
dominant  process  for  each  KRF  routine.  The  cost  of  the 
algorithms  in  floating-point  operations  that  compose  the 
TKRF  receiver  are  displayed  in  Table  1.  Note  that  for  a 
matrix  of  dimensions  /,  x  /2,  the  complexity  of  its  SVD 
computation  is  around  0(/1/2mtn(/1,/2))  [3]. 

Table. 1 :  TKRF’s  Computational  complexity  iiAfobating 


operations 


Condition 

Complexity 

KRF-1 

P>MS 

min(MDKJ  ,N)MDK]N  Ms 

KRF-2 

]>M1 

min (Md  K,  Ms  )Md  KMsM1 

KRF-3 

k>m2 

min  (Md,M1)MdM1M2 
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It  is  noteworthy  that  KRF-1  tends  to  be  more 
complex  than  KRF-2,  while  the  latter  is  likely  more 
complex  than  KRF-3.  One  of  the  reasons  is  that  there  is 
likely  reduction  of  the  input  data  volume  after  each 
Khatri-Rao  Factorization. 

IV.  PERFORMANCE  EVALUATION 

For  the  assessment  of  the  functionality  of  the 
proposed  receiver.  Bit  Error  Rate  (BER)  and  channel 
Normalized  Mean  Square  Error  (NMSE)  are  evaluated  . 

The  coding  matrices  C,  G^1^1  and  are  (tmncated) 

DFT  matrices,  and  we  assume  H(SRi;i  ~  CK (0,1/MS), 
H(«ir2)  ~  CK^.l/M^)  and  H(,?2D)  ~  CJV(0,1/M2). 
Symbol  energy  is  given  by  Es,  and  the  additive  noise 
samples  at  relays  and  destination  nodes  are  complex 
standard  normal  random  variables  wih  variance  equal  to 
one.  The  Bit  Error  Rate  (BER)  and  Channel  Normalized 
Mean  Square  Error  (NMSE)  curves  are  evaluated  using 
107  runs  of  Monte  Carlo  simulations. 

The  impact  of  the  code  length  P  and  J,  put  in 
practice  by  the  source  and  by  the  first  relay  station,  has 
already  been  investigated  in  two-hop  systems  ,  so  it  is 
adequate  to  assume  -  and  it  was  verified  -  that  the  KRST 
coding  at  R2  also  adds  a  new  coding  gain  on  the  outcome 
of  symbol  estimation  in  function  of  K .  Therefore,  perhaps 
it  is  more  interesting  to  evaluate  how  the  TKRF  receiver 
behaves  when  the  size  of  cooperative  network  increases. 

In  the  proposed  model,  it  means  when  Ms ,  M1,  M2  and 
Md  increase.  Fig.  3  demonstrates  the  BER  curves  for  M  = 

Ms  =  M1  =  M2  =  Md  =  {2,3,4}.  Ms  is  set  to  2  to  allow 
the  same  number  of  information  symbols  to  be  sent  in  all 
simulations.  Coding  spreading  lengths  obey  the  minimum 
values  to  satisfy  identifiability  condition  in  (15),  and  with 
N  =  50,  the  uniqueness  condition  in  Theorem  2  is  also 
complied.  For  two  modulations  (4-PSK  and  16-PSK),  the 
TKRF  receiver  responded  in  a  expected  manner  for  a 
MEMO  system:  more  antennas  meant  a  greater  diversity 
order  and  lower  BER  figures  as  Es  increase,  while  the 
higher  the  modulation  order,  the  higher  the  BER  value. 
According  to  Table  1,  the  estimation  complexity  from 
M  =  2  to  M  =  4  increased  almost  50  times.  This  huge 
increment  came  mostly  from  the  KRF-1  routine,  which 
corresponded  to  more  than  98%  of  the  overall  complexity. 
Whilst  this  number  does  seem  elevated,  there  is  still  a 
valuable  energy  gain  to  achieve  a  target  BER  (e.g.  «  10 
dB  at  BER  =  10“4  for  both  4-PSK  and  8-PSK). 
Moreover,  the  decision  to  let  the  receiving  node  do  all  the 
work  allows  the  use  of  the  AF  protocol  at  the  relays, 
saving  them  from  heavy  power-consuming  decoding 
processes.  Finally,  since  channel  estimation  (KRF-2  and 
KRF-3  routines)  offers  very  little  complexity  overall,  a 
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more  frequent  CSI  feedback  to  the  transmitting  nodes 
would  undoubtedly  allow  further  power  optimization. 

The  channel  NMSE  for  the  same  set  of 
simulation  parameters  is  shown  in  Fig.  4,  but  only  for  the 
4-PSK  simulation.  Besides,  the  curves  for  M  =  3  is 
omnitted  to  facilitate  the  visualization. 

P=MS«2,  M=J=K=Md=M1«M2 


The  slopes  of  the  curves  indicate  that  CSI 
recovery  responds  differently  for  KRF-2  and  KRF-3  in 
function  of  Es  -  KRF-1  is  responsible  for  symbol 
estimation.  For  H  <SRi\  approximately  an  order  of 
magnitude  of  the  error  is  reduce  per  decade,  while  for 
H(rir2)  and  H^2D^  the  KRF-3  algorithm  reduces  the 
NMSE  in  two  orders  of  magnitude  per  decade. 


V.  CONCLUSION 

This  article  brought  a  solution  for  the  issue  of 
blind  (symbol  and  channel)  estimation  in  one-way,  three- 
hop  non  regenerative  relaying  systems.  Based  on  -  the 
nesting  concept  of  the  4th-order  Nested  PARAFAC 
decomposition,  this  work  extends  this  idea  to  5th -order 
tensors,  allowing  the  leap  fromtwo-hop  systemto  a  three- 
hop  system  Alongside  the  proposition  of  the  new 
equations  for  this  new  derived  tensor  decomposition,  this 
work  brings  an  effective  non-iterative  semi-blind  receiver. 
Computational  simulations  validate  the  proposed  receiver 
and  bring  its  performance  in  terms  of  symbol  and  channel 
estimation,  as  well  as  complexity  and  transmission  rate. 
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It  is  important  to  regard  the  three-hop  system 
(and  its  receiver)  as  an  intermediate  step  between  the  two- 
hop  system  and  a  generalized  multi-hop  scenario.  There  is 
a  clear  pattern  that  indicates  that  several  expressions 
presented  in  this  paper  (as  the  identifiability  and 
uniqueness  conditions)  can  be  further  generalized  to  any 
number  of  hops.  However,  there  are  few  loose  ends  that 
need  to  be  addressed  in  the  next  works:  one  is  the 
complexity  one,  since  greater  number  of  hops  also  means 
a  much  larger  volume  of  data;  another  one,  now 
appreciated,  is  the  expressive  variety  of  possibilities  that 
opens  up  to  deal  with  the  data  and  to  solve  the  issue  of 
joint  symbol  and  channel  estimation  in  a  multi-hop 
scenario. 
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